VORTEX DYNAMICS IN THE PRESENCE OF EXCESS 
ENERGY FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION 

MATTHIAS KURZKE, CHRISTOF MELCHER, ROGER MOSER, AND DANIEL SPIRN 

Abstract. We study the Landau-Lifshitz-Gilbert equation for the dynamics 
On ■ oi a magnetic vortex system. We present a PDE-based method for proving 

vortex dynamics that does not rely on strong well-preparedness of the initial 

C ^ ' data and allows for instantaneous changes in the strength of the gyrovector 

Cn force due to bubbling events. The main tools are estimates of the Hodge 

J, , decomposition of the supercurrent and an analysis of the defect measure of 

C^ ' weak convergence of the stress energy tensor. Ginzburg-Landau equations 

with mixed dynamics in the presence of excess energy are also discussed. 
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1. Introduction 

Ferromagnets in a domain i7 C M'^ are typically modeled by a magnetization vec- 
tor m : O — > §^ with values on the unit 2-spliere. Dynamically, the magnetization 
satisfies the Landau-Lifshitz-Gilbert (LLG) equation: 

,^, dm ( dm 

(1) — =mxU— -/.eff 



where h^s is the effective field, arising from the L^ gradient of the micromagnetic 

^ , energy E^{m), and a > is a damping parameter. The form of E^{m) depends 

^O ' heavily on the physics of the ferromagnetic sample and contains a nonlocal pseu- 

CN , dodifferential operator; however, for thin, isotropic materials the energy simplifies. 

When the domain is thin, the magnetization vector lies mostly in the plane; con- 
sequently, 971 — (mi, ?Ti2, 013) =: (771,7713) behaves roughly like m k, (771,0) where 
^1^ ' 771 G S^. Furthermore, constraining m onto the plane induces the formation of 

vortices^ points (ai, . . . ,aAr) G fi^, about which the winding number of the pla- 
nar component, m^ is quantized. In the center of each vortex, the micromagnetic 
vector m k, (0,0, ±1); and hence, such vortices carry both the S^-degree of the 
winding number and a polarity, which may be interpreted as an §^-degree and is 
in simple situations given by the value of 7713 at the center of the vortex. Hence, 



JH ' micromagnetic vortices carry two pieces of information, whereas vortices arising in 

■ - - superconductivity, superfluids, and Bose-Einstein condensates have only the wind- 

ing number. 

The micromagnetic energy and the local area element concentrate at the site of 
the vortices (ai, . . . , ajv). In particular the local area element, described by the 
micromagnetic vorticity, concentrates via a cover of a hemisphere at the site of a 
vortex. The question of how such concentrated quantities behave in ([T]) has been a 
rich field of study. Thiele 13^ and Huber fTT! showed formally that in the dynamical 
setting concentrations obey an equation of the form 

Fn + Gn y- CLn- aQfln = 
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with an interaction force Fn — Fn{ai, . . . , un), a gyrovector force Gn = 4g„e3, and 
a damping term ao<in- 

In a previous paper j22l . the authors gave further justification of this motion law 
by deriving it rigorously from the LLG equation, but only under strong assumptions 
on the initial data and for boundary conditions, chosen mainly for technical rather 
than physical reasons. A similar problem with an additional spin transfer torque 
was studied in a different paper |20j by the first three authors. The undamped 
problem was studied by Lin-Shatah [3S] and later by Lin- Wei [35], who found 
traveling wave solutions for vortex pairs. 

Here we present a new approach that works under significantly weaker assump- 
tions and for other types of boundary conditions. While the previous work relied 
heavily on variational methods and ideas of Lin-Xin [37], Colliander-Jerrard [B], 
and Sandier-Serfaty [231 El], we now mostly use the structure of the PDE coming 
from (HI) and compensated compactness arguments in the spirit of Helein [16_ and 
Lin- Riviere [21] . 

1.1. Mathematical setting and results. We use essentially the same terminol- 
ogy and notation as in our previous paper [22J. Therefore, we keep their discussion 
brief. 

We use the functional 

(2) E^{m) = / e^{m)dx 

Jn 

with 

(3) e,{m) = ^\Vm\^ + ^ml 

as a model for the micromagnetic energy. Formally when e is very small and the 
energy is not too large, then m must take values close to the equator S^ x {0} in 
most of n. For topological reasons it may be forced to reach the poles at certain 
points, but the third component will decay rapidly away from these points. In 
particular one expects a sequence of local minimizers of ^ to converge to a map 
m^ — (m^, 0) away from a finite number of points ai, . . . , gn in Q. 

In the vicinity of each vortex the energy blows up at a rate tt log - up to lower 
order terms. Hang-Lin H^, making use of arguments of Bethuel-Brezis- Helein [2] 
for the associated Ginzburg-Landau energy, showed that as e — > 0, 

E,{m) =iv(7rlogi+7J +W{a) + o{l), 

where 7 is a universal constant and 

W{a) — — TT 2_. log |an — flml + boundary effects 

is a renormalized energy when the winding number about each vortex is one. The 
renormalized energy is the Kirchoff-Onsager functional arising in Euler point vor- 
tices. 

We allow winding numbers ±1 in this paper; therefore, we consider pairs (a„, d„) 
of points a„ G ri and winding numbers d„ = ±1 for n = 1, ... ,7V. Writing a = 
(ai, . . . , a„) and d = (di, . . . , d^), we obtain a renormalized energy W{a, d). This 
will also depend on the boundary conditions, and we study Dirichlet as well as 
Neumann boundary data. As the renormalized energy is the same as in the theory 
of Ginzburg-Landau vortices and is discussed extensively elsewhere [5], we do not 
write down the details here. 
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Due to the blowup of energy at the site of each vortex, one expects that 

I \ ^ 

log 7 ^1 

in distribution when we have a family of maps with good properties. Apart from 
the energy density, another fundamental quantity associated to the micromagnetic 
vector m is the vorticity, 

, , / dm dm\ 

which can be viewed as the Jacobian of the mapping m : J7 — > S^; consequently, 
it measures the local signed area of the mapping. If m is smooth on a disk Br (x) 
and '>TT'\gBr-(x) takes values on the equator S^ x {0} with winding number d, then 



L 



uj{m) dx £ 27rd + 47rZ. 

iBAx) 

Due to the "relaxed constraint" which induces m to lie mostly in §^, one finds 
that the vorticity will concentrate at a number of points, each point giving rise 
to a concentration of the vorticity of the amount 2'Kd + 47rq; namely the area of 
d hemispheres plus the area of a number of "bubbles" , or covers of §^ . Thus we 
expect 

N N P 

n—l n— 1 p—1 

where g„ G Z correspond to bubbles at the vortices and qp £ Z to bubbles elsewhere. 
Under strict conditions on the amount of energy, one can show that q„ = and 
qp = for all n and p, see [35] . 

In order to separate the winding numbers dn from the §^ degrees q„, we also 
consider the planar Jacobian 

dmi dm2 dmi dm2 
dxi dx2 0x2 dxi 
and we expect that 

N 
n=\ 

as well. Studying the Jacobian also helps for the analysis since we can apply 
standard results from the theory of Ginzburg-Landau vortices to Jim). 

The question of how concentrations, described above, are moved by the Landau- 
Lifshitz-Gilbert equation was considered by the authors in f35^. Using the model 
£*£ for the free energy and the abbreviation 

(4) f^{m^) ~ Am + |Vm| m (rn^e^ — rri^mj 

for the negative L^ gradient, the equation takes the form 

UTTi ( OTXX 

(5) -^ =me X U,-^-/^(m,) 

and the behavior of a^ for e — )■ is crucial for the answer. If it decays very slowly 
or not at all, then we expect no gyrovector force in the limit; on the other hand, 
if it decays rapidly, then the gyrovector force will dominate and the damping term 
a^an will be invisible. The most interesting case is when both terms coexist. This 
is expected when a^ is of the order 1/| loge|. We assume that 

ae log !> Qo as e — > 0. 

e 
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In the aforementioned paper [22, . a combination of differential identities for the 
time-evolution of uj{m^) and ee{m^), along with special choices of test functions, 
yielded an ODE for the evolution of the vortex positions: 

dW 
(6) 47rg„id„ + 7raod„ + - — (a, (i)=0, 

dan 

where g„ — ±|. The proof of ^ in [52] relies heavily on the so-called well- 
preparedness of the initial data. That is, we define the excess energy to be the 
difference between the actual energy and an expression that describes asymptoti- 
cally the energy needed to develop the observed vortices. It is then assumed that 
the excess energy tends to as e — >■ 0. Furthermore, the arguments require Dirichlet 
boundary data and the assumption that all winding numbers d„ = 2g„ are of the 
same sign. 

It is natural to ask what happens to concentrations in ^ when the initial data 
are not well-prepared, and there are two important problems that occur. The 
first issue is that we cannot prevent bubbling from occurring near vortices. In 
22] the authors prevented this by strict energy control, which no longer holds. 
Consequently, the gyrovector value can spontaneously change. The more troubling 
problem concerns the proof of the convergence of the stress energy tensor to the 
stress energy tensor of the limiting canonical harmonic map m* since energy bounds 
only imply weak convergence in H^ away from the vortices. To overcome the lack 
of strong convergence, the authors in [22] prove bounds for J^^^J^ ,^ -, iVm^ — 
Vm*p dx in terms of the excess energy; hence, if the initial data are well-prepared 
(i.e., the excess energy vanishes at the initial time), then m^ converges strongly to 
m* outside of the vortex cores. Again this argument fails when the initial data are 
not well-prepared. 

This paper develops techniques to establish the vortex motion law of micromag- 
netic vortices in the presence of excess energy. The arguments work when winding 
numbers of both signs are present and for both Dirichlet and Neumann boundary 
conditions. 

Let g = {g, 0) : dfl ^ S^ x {0} be a smooth map and consider smooth initial 
data m° : 17 ^> §^ with m^jao = g. We assume that there exists a constant Cq 
such that 

E,{m°) <7TNlog- + Co 
e 

for every e e (0, |]. Furthermore, we assume that there exist a° = (aj, . . . , a%) G 
U^ with a„i 7^ a-ii for m ^ n, d = {di, . . . , dN) with dn = ±1, and g" — {qi, . . . , q%) 
with (/J^ e 4 + Z, such that 



log 7 



TV JV JV 

-^ TT y^ SaO , J{jnl) ~^ ^ / dnSaO , and a;(m°) -^ 47r \^ Qn^a^ 



in distribution. 

These are the assumptions we use for the problem with Dirichlet boundary con- 
ditions. When we study Neumann data, then we drop the condition m^lgji — g, 
and all of the other assumptions remain the same. 

Theorem 1.1. For every e G (0, 2], there exists a weak solution m^ of ([5]), with 
mg(0) = mg and m(t)\Q^ — g for all t > Q, that is smooth up to finitely many 
points in space-time. Furthermore, there exist a number T > 0, a piecewise constant 
function q : (0,T) — >■ (^ + Z)^ with finitely many jumps, and a sequence efe — )■ 
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such that for almost every t £ (0,T), 

N N 



ee, ("!•.(*)) 



and 






N N 



n— 1 n— 1 



where a — (oi, . . . , ajv) solves 



dW 
Airqnidn + TraoOn + t; — (a, rf) = 0, n = 1, . . . , A^, 

TO (0,T) w;ii/ia(0) = 0°. 

T/ie same statement holds for homogeneous Neumann instead of Dirichlet bound- 
ary conditions. 

One unusual feature of the vortex motion law is the possibility of spontaneous 
changes in the gyrovector, represented by the jumps of g„. For a single vortex in 
a disk this would correspond to a sudden change in the direction or the speed of 
circulation of a vortex about the origin. 

A consequence is that using only Theorem ll.il we cannot predict the trajectories 
from the initial data alone. This is also the reason why we make a statement only 
for a sequence e^ — > 0. Unless it is possible to extract more information, it is 
conceivable that for a different choice of e^, the trajectories of the vortices are 
different, and in such a situation we would not have convergence of the entire 
family m^. 

From a physical point of view, changes in the direction of rotation are not some- 
thing completely unexpected. They are known to occur in field- or current-driven 
switching processes, where a vortex is driven from its equilibrium position by an 
applied field or current (spiralling outward), then changes its polarization and spi- 
rals back into equilibrium after the external field or current has been switched off, 
see e.g. [TTl [T3] for some related simulations. 

In order to prove the theorem, we need to establish two things. First we show 
that the stress energy tensor Vm^ (g) Vrrie converges to Vm* (g) Vm* outside of 
a defect measure concentrated at a set of delta functions. Second we analyze how 
the defect measure interacts with the micromagnetic vortices by applying different 
test functions on the differential identities for uj{m^) and ee{m^), and it is shown 
that the defect measure has no effect on the vortex dynamics. 

1.2. Ginzburg-Landau equation with mixed dynamics. The Landau-Lifshitz- 
Gilbert equation given above describes a mostly-planar §^ micromagnetic vector. 
There is an analogous problem in complex Ginzburg-Landau theory in which an 
order parameter u : M" — > C evolves according to an equation with mixed dynamics: 

(7) ^a, + ^)^ = Au+l^u{l-\u\'), 

which is also a hybrid of gradient flow and Schrodinger dynamics. Smooth solutions 
to ([7]) satisfy an energy dissipation equality 

2 



E,i{u{t)) 



du 
'dt 



dxds = £^gi('u(0)). 



Jn 

2 , 1 



where Eg\{u) = J^ Gg\{u) dx and egi(u) = ^ | Vu| +4^ (l — |up) • In the Ginzburg- 

log - 



Landau setting, a vortex is defined by the concentration of energy "f" " — )> ^ 7r5a„ 



6 M. KURZKE, C. MELCHER, R. MOSER, AND D. SPIRN 

and the concentration of the Jacobian J{u) — >■ "^ndnSa^, both of which are well 
understood, see [TS] , 

The formal vortex motion law of E [8] was established rigorously by Miot [28] in 
the plane and the authors of this paper |21j on bounded domains. In both papers the 
authors again use the strong well-preparedness of the initial data to show that the 
stress-energy tensor converges to the stress-energy tensor of the limiting canonical 
harmonic map as e — >■ 0. This trick, pioneered by Colliander-Jerrard [6j and Lin- 
Xin ^7\ in the context of Ginzburg-Landau vortex dynamics, is reliant on excess 
energy estimates. One drawback of this method is that the initial data must lie 
close to the optimal mapping in an H^ sense. This certainly fails to take advantage 
of any energy dissipation that is present in ([7]) , as is well understood in the purely 
dissipative problem. 

When the dynamics are purely dissipative, much stronger results can be shown 
that remove the well-preparedness assumption on the initial data. The first proofs 
of the vortex motion law for the Ginzburg-Landau heat equation a^ ^&^ = Aue + 
■i^Ue (l — |uep) by Lin [53] and Jerrard-Soner ^19] proved strong convergence of 
Vue -^ Vu*, the limiting canonical harmonic map, away from the vortex cores by 
parabolic estimates, even when the excess energy is of order 0(1). The proofs de- 
pend on using estimates of the form Au^-l-^Ue (l — \u^\^) = 0^(1), which arise from 
energy bounds. Later refinements of the vortex motion law by Bethuel-Orlandi- 
Smets [5] and Serfaty |34j again used the purely dissipative nature of the equation. 
A naive adaptation of these methods fail for ([7]) : while we have the same L^ control 
of ^^ as in the purely dissipative case, it is not possible to deduce smallness of 
Auf + -^Uf (l — |uep). The reason for the difference is that jag + i| — > 1 while 
Ickel ^- as e -^ 0. 

For initial data without vortices and close to a constant map, Miot [3^1 studied 
([7)1 in the damped wave regime. It is unclear how this approach can be generalized 
to initial data with vortices. 

Using similar arguments to the proof of Theorem ll.il we can establish the vortex 
motion law for ([7]) for initial data that are not well-prepared. 

Theorem 1.2. Let Ue{t) be a sequence of solutions to ([7]) with initial data Ue{0) = 
Wg and either Dirichlet or homogeneous Neumann boundary data. Assume that 
Ee{u'^) < A^Trlog -^ + Co for some constant Co and 

eJu°) A ^ 

— — -, > y T^SaO and J{u^) — ^ / T^dnS^o 

log - ^-^ " ^-^ 

° t n=l n=l 

with dn — ±1 for n — 1, . . . , A^. Then there exists T > such that for all t £ (0, T) 
( (t)) ^ ^ 

(8) " \ ^^^^a„{t) and J{u^{t)) -^ ^TTdnSa„(t), 

° e n=l n=l 

where a„(i) solves 

dW 

(9) 2'Kdnian + TTaoCLn + -^ — (a, d) — Q 



dOn 



for n = 1, . . . , N . 



Therefore, the vortex motion law is exactly as in |21| . and, unlike in the micro- 
magnetic case, there are no changes in the direction of rotation of a single vortex 
about the origin in disk domains. Since we are unable to handle vortex collisions or 
a vortex migrating to the boundary, the T in the theorem represents the first time 
for which either of these two events occur. Finally, we note that the limiting motion 
law ^ is independent of subsequence and uniquely determined by the initial data, 
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unlike in the LLG case. The proof of Theorem 11.21 is very similar to the proof of 
Theorem 11.11 and so we only briefly discuss the proof in section [T] 

2. Mathematical tools 

Here we provide some mathematical notation, discuss topological quantities aris- 
ing in our study, and review facts about the expansion of the micromagnetic energy. 

2.1. Notation. We first introduce some notation that helps keep track of the po- 
sitions of the expected vortex centers. We typically write a = (ai, . . . , qn) £ il^ 
for these positions. To control collision and escape to the boundary, we define a 
minimal measure of intervortex and vortex-boundary distance: 

p{a) = min < — min ja™ — a„| , min dist(a„,9r2) 

I 2 7n^n n—l,...^N 

As we expect the energy to concentrate near the vortices, we sometimes need to 
cut out vortex balls of radius r or merely the vortex centers. We define 

N 

ftr{a) := il\ [J Brian), ^o(a) := ^\{ai, . . . , un}- 

We denote (•, •) to be the scalar product for elements in R^, and let (•, •) denote 
the R^ scalar product, which will arise often when we project R'^ onto R^ x {0}. 

2.2. Energy density, vorticity, and the planar Jacobian. The analogous 
quantity of uj(ra) for u : R^ ^' C is the Jacobian J(u). The Jacobian is useful 
for describing the winding number about a vortex in Ginzburg-Landau theory. We 
can also write J{u) = ^ curlj(w), where the supercurrent j{u) is defined as 

j(u) = {iu,Vu), 

for the complex-inner product (•, •). The connections between uj{m) and J{m), for 
the projection of m onto R^ (and complexified) can be seen by: 

J{m) — m,zuj{m) 

uj(rn) — 3?7i3 J(?7i) + curl {m2'm^\Imi — ?nim3VTO2) , 

both of which follow from direct calculations [20l Lemma 1] . 

Our method of proof entails tracking how concentrations in the energy density 
and in the vorticity move in time. In that regard we write down differential identities 
that are satisfied by ee(m) and uj{m). In particular we have: 

d 2 

(10) —ee{m) = a^ {f^{m),\/m) - div (m x f^{m), Vm) - a^ |/e(m) 

and 

d 

(11) — w(m) = curl(/g(m), Vm) + a^ curl (m x f^(m), Vm) , 

where f^ is defined in (H]). See ^2; for derivations of these differential identities. 

2.3. Renormalized energy. The micromagnetic energy E^(rn) can be expanded 
up to second order. This expansion, presented in Hang-Lin [13], says that the 
micromagnetic energy behaves roughly like 

E,{m) K. N{t: log -+-/) + W{a, d), 
e 

where W{a, d) is the renormalized energy of Bethucl-Brezis-Helein [5]. 
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To define W{a, d) we first write down the canonical harmonic map m^ : fJ — > §^, 
which is defined via the Hodge system 

N 

divj(m*) = and curl j("^*) — / ^ 27r(5a^ 

n=l 

with cither the Dirichlet condition 

m, = g 
on dfl or the Neumann condition 

on dVl, depending on the problem that we study. Then m* can be written as 



m* = e 



N 
iS TT I X ~ Or, 



n 



L via; - a„| 



for a harmonic function : J7 — > M satisfying the appropriate boundary conditions. 
Then 



W{a, d) = lim 



/ ^|Vm,|^-iV7rlogi 

= — TT 2^ djndn log \am ~ On | + boundary effects, 

which is precisely the Coulombic energy of point particles (with a repulsive bound- 
ary in the case of a Dirichlet boundary condition). 

Finally, we note that if e C^{fl) such that V^V0 vanishes near a vortex, 
then we have [5] 

^ d r 

(12) TT V V^(/)(a„) • T—W{a, d) = - / V-^Vcf) : (Vm, ® Vm*) dx. 



This identity will be important for the derivation of the motion law of the vortices. 

3. The Landau-Lifshitz-Gilbert equation 

The proof of Theorem 11.11 requires a combination of ideas from various theories 
and spans several sections. Most of the arguments do not make use of the boundary 
conditions at all; therefore, we treat Dirichlet and Neumann data simultaneously. 

In this section, we discuss properties of the Landau-Lifshitz-Gilbert equation 
that can be described as well-known, even though they may be difficult to find in 
the literature for the exact version of the equation used here. These properties have 
been derived for related equations — in particular the harmonic map heat flow [551151 
[71 [Sni a-iid a simpler version of the Landau-Lifshitz-Gilbert equation '12] . Examining 
the corresponding arguments, we see that they carry over to our situation. 

We fix e > for the moment. Suppose that we have initial data mo G C°°{n: §^), 
and we want to solve the Landau-Lifshitz-Gilbert equation under either Dirichlet 
or Neumann boundary conditions. Following Struwe [35^, Chang [5], and Guo- 
Hong [12], we can construct a weak solution m with m(0) = mo that is smooth 
away from finitely many points. That is, there exists a finite set of singular points 
S = {(ti,a;i),...,(i^a;0} C (0, cxd) x 17 such that m e C°°(([0,oo) x f7)\I];§2). 
Furthermore, the map m belongs to the Sobolev space 

i/i((0,r);L2(f7;S2))nL°°((0,oo);i/i(f7;S2)) 

for every T > 0. 
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The singularities arising in this construction are a consequence of energy con- 
centration. If (i*,a;*) is a point in S and r > is small enough to separate this 
point from all other singularities, then we have a g^ G i^% such that 



/ e(:(m,) dx + 47r |qj| < liminf / ee{m)dx, 

uj{m) dx + 4:TTqi = lim / uj{m)dx. 

X') i^t^ J{t}y.Br(x^) 



' {t^]y.Br(x 

Indeed, the number qi has a geometric interpretation, and in order to understand 
it, we need to analyze the process leading to singularities. This is done with the 
tools developed by Ding-Tian 7 and Qing [30 . 

Fixing the singular point (t*,x*), we can find sequences Xj — t- x*, tj /^ t\ and 
Tj \ 0, such that the rescaled maps mj{x) = m(tj, tjx + Xj) converge to a critical 
point of the Dirichlet functional. These critical points are called harmonic maps, 
and if they arise through a blow-up process as here, we also speak of harmonic 
bubbles. Via the stereographic projection, every harmonic bubble at a singularity 
in the interior of J7 is identified with a harmonic map S^ — > S'^. Since such harmonic 
maps are very well understood ||lj, we have good information about their energies. 
If q is the topological degree of a bubble at a singularity in the interior of il, then the 
corresponding energy is 47r|q'|. Under Dirichlet boundary conditions, this is the only 
situation possible. Under Neumann conditions, bubbles at the boundary are also 
conceivable. They give rise to harmonic maps on a hemisphere (with Neumann 
boundary conditions on the equator), which can be extended to harmonic maps 
on §^ by reflection. Thus we can think of them as half-bubbles, and the degree 
is naturally measured by a half-integer q. With this convention, the energy of a 
bubble remains 47r|g|. 

The number q^ in the preceding formulas is the sum of the degrees of one or 
several bubbles at the singularity (t*, x^). In the first formula, we have an inequality 
only, because degrees can cancel each other, whereas the energy is always positive. 
Equality holds if all the degrees have the same sign. In any case, we have at least 
one bubble and therefore 

/ Ce (m) dx + 27r < lim inf / e^{m)dx. 

J{t'}xB,.{x') */'*' J{t}xB,.{x') 

Using the formula at every singular time and the energy conservation law (|10p 
between singular times, we obtain the energy inequality 

dm 



(13) E,im{t2)) + 2TT\Y.n{{ti,t2] xn)\+a, 



U r- a___ 2 



dt 



dxdt < E^{m{ti)) 



for < ti < t2. 

Weak solutions of the Landau-Lifshitz-Gilbert equation are not necessarily unique 
[II 137). but uniqueness follows when we impose suitable additional conditions. In 
particular, it follows from arguments of Freire [10] (see also Harpes J^ and Rupflin 
[3T]) that there is exactly one weak solution such that the energy is non- increasing 
in time. We call it the energy decreasing solution. 

We now consider initial data m° S C°°(S1;§^) that satisfy the assumptions in 
section 11.11 We study the energy decreasing solutions mg of the Landau-Lifshitz- 
Gilbert equation ([5]) in (0, oo) x O with 171^(0) = Tn°. Then each m^ is smooth away 
from finitely many singular times. That is, there exists a finite set S^ C (0, cxd) such 
that m^ is smooth in ([0, oo)\I]£) x il. Furthermore, we have the energy inequality 

2 



S,(m,(i2)) + 2^|E, n(ti,t2]| + a, / [ 

Jt, Jn 



dm. 



dt 



dxdt < E^{m{ti)) 
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for all ti,t2 G [0,00) with ii < ^2- Thus we have derived the first statement of 
Theorem 11.11 from known arguments. It is more difficult to prove the motion law 
for the vortices, although in the first step, we still rely mostly on known results and 
ideas. 

4. Continuous vortex motion 

The aim of this section is to show that the vortices (centered at a^, . . . , a^ at time 
0) persist for positive times and the vortex centers have a continuous trajectory. 
More precisely, we want to prove the following. 

Theorem 4.1. There exist a sequence e^ \ 0, a number T > 0, and a continuous 
curve a : [0,T) — )> il^ with a(0) = a^ , such that p{a(t)) > and 

N N 

ae^e,, (m,, (t)) ^ tt ^ 5a„(t) and J{m,^ (t)) ^ tt ^ c?n'5a„(t) 

n— 1 71=1 

as distributions for every t e [0,T) as k -^ cxd. 

Moreover, there exists a finite set S C [0, T] with the following property. Suppose 
that [^1,^2] C [0,r] is an interval with En[ii,i2] = 0- Then there is a number K G N 
such that m^^ is smooth in (ii,t2) x Q for every k > K. 

Proof. The convergence of the Jacobians Jim^^ ) (for a suitable subsequence) follows 
almost directly from results of Sandier-Serfaty |33j . Namely, Theorem 3 implies that 
there exist a sequence efe \ and a function J„ G C°^^/^([0,T); A^(f7)) such that 

J(m,,(i))^ J,(i) 

in W~^'^{il.) for every t G [0,T). On the other hand, the logarithmic energy bound 
obtained from the energy inequality p^ . together with results of Jerrard-Soner 
[TBI Theorem 3.1], give good information about the structure of J*(i): it must be 
of the form 

N 

J,(i) = 7r^rf„(t)4„(t) 

n=l 

for certain integers di(t), . . . , dff{t) and certain points ai(t), . . . , a^(i) G fl. Let 
77 > 0. By the continuity of J,, we can choose T small enough so that 



N N 

n— 1 71—1 



w-i-i-{n) 



If rj is chosen sufficiently small, then it follows that N ~ N and dn{t) = dn for 
n = 1, . . . ,N. Moreover, the continuity of J^, also implies that the curves oi, . . . , on 
are continuous in [0,T). 

The convergence of the Jacobians implies that for every t G [0, T) and for every 
r G (0,p(a(i))], we have 

/ Ce Jme J da; > TT log C. 



Hence 



/ ee^.(mejrfa; -> 0, 



'Sl^(a(t)) 

and the convergence of the energy density follows as well. 
For every t G [0,T), we have 

liminf ( i?£^(mc^(t)) - A^TT log — ) > -00. 
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We may assume that the same holds at t = T (otherwise we replace T by a smaller 
number). Then we have a uniform bound for the number of singular times in [0, T] 
by (fT3l) . and we may assume that there exist a finite set E C [0,r] such that 

limsup max min Is — tl = 0. 

fc_i.oo sGS.j^n[o,T] tes 

This set then has the required property. D 

5. Strong convergence 

In this section, we study the behavior of the solutions at a fixed time where no 
concentration of the kinetic energy occurs. It is then convenient to rescale the time 
axis, and therefore we now consider the family of equations 



We assume that we have weak solutions m^ G C°°([0, 1) x 17; §^) satisfying 

(15) lim / / ^ dxdt^O. 

Let a € n^ and suppose that there exists a function 7 : (0, 00) — > M. such that for 
any r > 0, 



limsup sup / ee{m^{t,x)) dx < "f{r). 

e-*0 0<t<lJn^{a) 

Furthermore, we assume that m^ -^ m* weakly in Hl^^{{0, 1) x ilo(a)) and strongly 

Lemma 5.1. Under these assumptions the following holds true. 

(i) The limit is a time-independent harmonic map into S^ x {0}, i.e., it sat- 
isfies 771,3 = and !^* — 0, and 

divj(m,) — 0, curlj(r7i,) = in (0, 1) x flo{a). 

(ii) Consider the measures 

fi.^C'^l. (Vme(O) ® Vm,(0)). 

Let r > 0. There exist a sequence efe \, 0, finitely many points bi, . . . ,bp G 
fir {a), and matrices -Bi, . . . , Bp G M^^^ such that 

p 
fi,^ ^ C^i- (Vm, ® Vm,) + Y^ BpSb^ 

p=i 

weakly* in (Cq (ilr(a)))*. Furthermore, the number P depends only on 

7(r). 

Remark. We will apply this result in a situation where the assumptions are justified 
only for a subsequence e^ \ 0. We give a proof only for the result as stated in 
order to avoid awkward notation; however, everything works for a subsequence as 
well. 

Proof. It follows immediately from the assumptions that to,3 = and ^^^ = 0. 
As a map with values in a 1-dimensional manifold, the limit automatically satisfies 

curl j(m*) = in (0, 1) x floia). 

Equation (fH)) implies 



1 (JTITi "^ ( (JTfl 

(16) divj(?Tie) = — ^ — V im^, 



dt \ " dt 
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Multiplying by a test function and integrating by parts in the first term on the 
right hand side, we also find 

div j(m*) = in (0, f) x fio(a)- 

ft follows by elliptic regularity theory that m» is smooth away from the vortices. 
Let S > and suppose that we have a disk B2r{xo) C fl with 



(17) 



limsup / ee(me(0))(ix < S^. 



We first want to show that a similar inequality persists for some positive time. To 
this end, let C e C^{B2r{xQ)) with < C < 1 and |VC| < C/r. Here and in the 
following, we indiscriminately use the symbol C for various universal constants. Set 

2 

dx. 
Then we have 

2 



F,{t)= f C^e,{m,)dx and G,{t) = f 

J{t}xn J{t}y 



dm. 



dt 



Kit) = - f 
J{t} 



e 



xO 



dm. 



dt 



C 



<-VFJt)GS)- 
r 



It follows that 



Hence 



dx~2 I CVC 

'{t}xO 



C 



drrit 



, Vrrif ) dx 



-VW)<-VGS)- 



^/K{s)<c5+-(s r [ 

"f" y Jo Jn 



dm. 



dt 



1/2 



dx dt 



and in particular 



f ee{m^{t,x))dx<CS^ 

JBr-ixo) 



r(a:o) 

forie [0,r2j2]. Set T = r2 ,52. 

Choose a cut-off function r; G C^{Br{xo)) with < r/ < 1 and IV77I < C/r. 
Consider the Hodge decomposition 



V{j{m,{t)) - j(m,)) = V(/.e + V^^e 



then 



|lV(/.e(t)|li2(R2) + |lVV^,(i)Ili2(R2) = h(jK(t)) - j("^*))lli2(B„(.o)) 

for every t e [0,t]. Using ((T6)) . we also obtain 

1 
^e Jo Jn 



Jn 



ri^\j{m^)~ j{m^)Y dxdt^ ^ I / 770, f dx dt 

a, . . o dt 



Jn 



, / dm^\ 

ricpe im^, —- — dxdt 
Jn V dt J 

<t>e^'n' (j("^e) ~ j{'m.*))dxdt 

r^V^Ve ■ UiiTT-f-) — Ji'm^,))dxdt. 
Jn 

For almost every t £ (0,t), we have (j>e{t) — > strongly in L^^^{M?) for every 
p < 00. Using Lebesgue's convergence theorem, we prove that 4',i{t) -^ strongly 
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inLf^^((0,T) x]R2)_ Hence 



and 



Jn 



(f'e'Vr/ ■ {j{m^) — j{m.f))dxdt -^ 



Jn 



rj(j)f^ im 



dm,/: 
~~dt 



dxdt ^0 



as e ^ 0. Moreover, 

ryV Ve • (Ji^TT-e) — Ji^TT-*)) dx dt 



Jn 



Jn 



r;V ?Ae • (zTTie — ITOh., VTOc) dx dt 



+ / V ipe- {i'm*,y{rj{mf_ — ■m^)))dxdt 

Jo Jn 



V ip ■ \/ri{im^,m^ — m^) dxdt. 



Jn 



We clearly have 



'0 Jn 
For almost every t, we observe that 



/ \7 ip ■ \/r]{im^,,m^ — m^,)dxdt ^ 0. 

Jo Jn 



C. 



[riiim^ - iTO*)]BMO(K2) < C||ry(VTOe - '^rn^)\\L^n) + —\\me - m*||i2(B^(^^)). 

The last term on the right hand side converges to 0. With compensated compactness 
arguments [3], we now estimate 

/ / 7/V ip^ ■ {im^ — irriif^'^m^) dxdt < C5 I i] iVm,: — 'Vm^,\ dxdt + o^{l). 

Jo Jn Jo Jn 

[im*]BMO(B,(2-o)) < CS, 



Similarly, 
and therefore. 



/ / V^'p^-{im^,'\/{ri{m^-m^)))dxdt <C5 I / yy^lVme- Vto*|^ da; di+Oe(l). 
Jo Jn Jo Jn 

Finally, we compute 



rjcf)^ " dxdt 
Jn dt 



riip^m^^ dx — I ri(j)^m^^ dx 

{t}xS1 J{o}xr2 



Jn 



rj—^m^sdxdt. 
at 



The first two terms on the right hand side are of order 0(e). Now we use the 

identity 

9 ., , „ f drrie ^ \ „ ( dm^ 
—j(m^) = 2 I— — ,Vme + V im^, ■ 



dt' 



dt 



dt 



hence, 
A 
We conclude that 



drrif 

dt n^"''^ 



, , dm^ T-^ \ ^ f 9m f 
div ( 2ry ( i——,\/m^ I - V77 I im^, — - 



dcp. 



dt 



<C 



LP((0,r)xR2) 
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for every p < 2. Thus 



1 



"e Jo Jn 
Combining these estimates, we obtain 



rjSg — - — dxdt -^ 0. 
' dt 



/ j Tf\j{m^)- j{m^)\'^ dxdt<C5 j JTf\\/m^~\Im^\'^dxdt + o^{l). 
Jq Jn Jo Jn 



Now let 



A, = <^ i e (0, r) : 



dm. 



dt 



it) 



< 



L2(a) 



Fix e G (0, ^] and t E A^. For a given xi E Br{xo), consider 



Then we have 
and 



rfi{x) = Tn^{t, ex + xi). 
||/i(rn)||L2(Si(o)) <Cae 



l|ei(m)|Ui(Bi(o)) <C6\ 

Choose /3 > 0. li 6 and e are sufficiently small, then known estimates from the 
theory of harmonic maps (e.g. [71 Lemma 2.1]) imply m^ < (5. We conclude that 
TO^3 < /3 in Ae X Brixo). 
Using the equation 

Ame3 = -|VmepTO£3 + — (me3 - m^g) + fe?,{m^), 



we obtain 



/ rj^\SIm,-i\'^dx^ j ifml^\Vm,\^dx-\ I 

J{t}xn J{t}xn £ J{t} 

-/ ri'^me3f^3{m^)dx -2 i]m. 

J{t}yA1 J{t}yA1 



rj (m,3 - m.a) da; 



V?7 • V?7ig3 dx. 



For i G Ae, we have 



/ rfml^\\/m^\'^ dx <2 j r]'^ml^\Vra^ ~Vm^\'^ dx 

J{t}y.n J{t}xn 

+ 2 / ?7^m^3|Vm,pda; 

< 2/3 / Ty^lVm^ - Vm*pdx + o,(l) 

i{t}xr2 



and rn^g — rn^g > |?7i^g in {t} x _Br(a;o). Furthermore, 



/ Tfm^dfez{m^) dx < Ca^ 

J{t}xn 



and 



Note also that 



1 



-2 / -qm^'i'S/ri ■ Vm^g dx < Celog 

'{t}xn ^ 



/ / rf'e^{m^)dxdt < Ce^ ( log- 

/" f v^ (\Vm^3\^ + ^] dxdt<Cf3 [ f rf\Vm,-\/m4^dxdt + o,{l). 
Jo Jn \ £" / Jo Jn 



It follows that 
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Combined with the previous estimate, this yields 

I if hj(m^)-j(m,)P + |VTOe3p + ^) dxdt 

<C{I3 + S) f f 7/2|Vm,-VmH.p dxdt + 0,(1). 
Jo Jn 

Since 



TTh 

Vtoc - Vm, = -j — —Vdl ~ m^3 + -j '^—{j{mf) - j{m^,)) 



■i ( 1 — 17 -m^ ) ii'm*), 



we have 



/ / ?7^|Vme — Vtti, p dx 



<C / ry2(|j(m,)-j(™*)l' + |Vm,3n dxdt + o,(l). 

Jo Jn 

For the integral over (0, t)\A^, we have the same kind of estimate as before. Thus 
/ rj'^ (\Vm,~Vm^f + ^j dxdt 

<C{I3 + S) [ f 7f\\/m, ~ \/m^\^ dx dt + o,{l). 
Jo Jn 

If /3 and S are chosen sufficiently small, then it follows that Vm^ — > Vm* and 

m^a/e -^ strongly in L^((0,t) x Br/2{xo))- 

Finally, let x G C^(Br/2{xo)) and consider again the energy identity from (jlOp 



■ / x^e,(me)dx ^ - / 

J{t}xO J{t} 



dm. 



dt 



2 



dx 



2/ xVx-(^i7^,Vm, ) dx. 

{t}xa 



Integrating in time and using (jlSp . we see that 

: / X^e.(.n.)dx^O as . ^ 



osc 

0<t<T 



Hence we have 



limsup / x^ec(me(0))dx < - / x^lVm^^dx. 

e->0 Jo 2 Jf2 



/o ^ Jn 

That is, we have m^{Q) -^ m* strongly in H^{Br/i{xoj) and me3(0)/e ^^ strongly 
in L'^{Br/4{xo)) — assuming the smallness condition (ITTI) . 

The second statement of the lemma now follows by a standard covering argument. 

D 



6. The motion law 

We consider the solutions m^ of ([S]) constructed in section |3] again. 

Let T > be the number from Theorem 14.11 We consider the sequence e^ \ 
from the same theorem and we use the same notation again. In particular E C [0, T] 
is the accumulation set of all the singular times of m^^ for fc ^^ oo. As E is finite 
and the vortex trajectories are continuous, it suffices to verify the motion law away 
from E. Let m* = (m*,0) : (0,T) x fi — > EP be the map such that m^,{t) is the 
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canonical harmonic map with vortices of degree d„ at a„(i), n ~ 1, . . . , N. Note 
that we may assume rriej. -^ m, locahy in H^ away from the vortices by part Q 
of Lemma [5.11 f apphed to time-rescaled versions of me). 

Consider a time interval (^1,^2) with E n [^1,^2] = 0, and so short that for some 
r e (0,/9(a(ii))], we have a„(t) £ B,^/2{an{ti)) for all t e [^1,^2] and n = 1, . . . ,N. 
Define the measures 

fle= jO^\- (Vme Vm,) 
on [ti,t2] X f2r(a(ii)). There exists a subsequence such that we have the convergence 
/le^ ^' fl for a matrix valued Radon measure /2 on [ii,i2] x ^r(fl(^i))- For almost 
every t G (^1,^2), there exist a measure /^(t) on flr{a{ti)) such that for all rj G 

(18) /" rjd^iit)^\hnj- f r^dfl. 

Jnr{a{ti)) ''N"" ^" J{t-h,t+h)xn,^(a{ti)) 

Moreover, for all S, G Co([ii,t2] x f2r(a(ii))), we have 

^dp,= ^ dfi{t) dt. 

I[ti,t2]xflr{a{ti)) Jti Jn^iaiti)) 

We also consider the measures 

H^{t) = £2 L (Vm^i) (K) Vm,(t)). 

Lemma 6.1. There exist finitely many functions bi,...,hp : [^1,^2] — ^ ^r(a(^i)) 
anrf Bi, . . . , Bp : [ti, ^2] — ?> M^^^ smc/i that for almost every t G [^1,^2], 

p 

fi{t) = ^^t) +^Bp{t)Si,^(^t)- 
p=i 

Proof. Fix t G (^1,^2) such that ^TE\i holds true and in addition, 



/ 



Crt+h 



dm. 



dt 



dx dt 



0. 



We may assume without loss of generality that the latter is true for almost every 
t, for a subsequence with this property exists. 

Let (?7£)£gN be a dense sequence in CQ(J7,.(a(fi))). Choose a sequence hk \ 0, 
such that 



1 
2^. 



k J[t-/ifc,t+/ifc]xnr(a(ti)) 



id^i 



m d^i) 



a,(a(ti)) 



< 



for ^ = 1, . . . , fc. We can now replace {mt^)k£fi by a subsequence (but still use the 
same notation), such that 



I j-t+hk 

2^/c Jt-hk 



Vi^tn^ 



Vmej. dx dt — 



for i = 1, . . . , fc. Finally, by the continuity of the functions 



Vt dn{t) 



< 



t^ r]eVm^^{t)<S>Vm^^{t)dx, 

Jn,.{a(ti)) 

we can find tk ^ t such that 

/ r]e\/m^^{tk)<S>Vm,_^{tk)dx- rjed^it) 



< 



for e^i,. 



,k. 



-a? 
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Define 

rhk{t,x) = m^jaeji- tk),x). 

Then rhk solves (fH|) for e^ instead of e, and tfie sequence satisfies the counterpart to 
condition (jlSp . Using Lemma lS.f I and passing to a subsequence again if necessary, 
we conclude that there exist 6i(i), . . . , fep(i) £ r2r(a(ti)) and Bi{t), . . . ,Bp{t) G 
]R2x2^ such that 

/ r]idfj.{t)^ mdfJ.*{t)+'^rie{bp{t))Bp{t) 

for all € G N. Furthermore, the number P is bounded by a constant independent of 
t by the energy estimates away from the vortices. The claim now follows. D 

Let (/>,V' e C^{^). Then we compute, using ^ and ((TTl), that 

J {t}xn 

- [ {alV^ + V^4>) ■ (f^{m,),Vm,) dx 
J{t}xn 

J{t}xn 
in [^1,^2]- Note also that 

/ V^(/)- (/,(m,),Vm£) dx= - / V^V0 : Vm^ ® Vm, da; 

J{t}xn J{t}xn 

by an integration by parts. 

First consider the identity with ■0 = and (p e C^{ilr{a{ti))). Define the 
measures 

A,(i) ^C'^\_uj{m^). 

Then the functions 

t^ / (/)(iAe(i) 

are uniformly Lipschitz continuous. It follows that there exists a sequence efc \ 

such that Ae^(i) ^ A(i) weakly* in (Co(r2r(a(ii))))* for almost every i, where X{t) 

is a Radon measure on flr{a{ti)). Furthermore, we can choose this subsequence 

such that the previous statements hold true. Lemma |6 . 1 1 then implies that X(t) is 

of the form 

p 

p=i 

for certain functions cti, . . . , ap : [<i, i2] ~^ K- Moreover, using Lemma 3 in [22| . we 
see that (7p{t) G 47rZ for every p and almost every t. 

Still using (J19p with -0 = and G C^(rir(a(^i))), we now obtain 



p=l 



p=i 
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Moreover, 

ti Jn 

since m^,{t) is a harmonic map. It follows that 

p 

p=i 

belongs to W^^°°{ti,t2) (so in particular it is continuous) and for almost every i, 

(20) J^ E ^vmiMi)) = E V^V</>(6p(t)) : B,{t). 

p— 1 p— 1 

We want to show that \{t) is constant, inserting test functions (j) with V^(/)(&p(i)) = 
0. However, since the Upit) may have a sign, the continuity of the above function 
does not immediately rule out nucleation or annihilation of dipoles. So we have to 
be a bit careful here. 

Define #(i) = |suppA(t)|. Choose io G (^1,^2) with #(to) = maxtj<t<t2 #(i) 
(which exists because the function takes integer values). 

Lemma 6.2. There exists an open interval I with to e /, such that #(t) = #(io) 
in I . 

Proof. Choose p > such that the balls B2p{x) for x G supp A(io) are contained in 
17 and pairwise disjoint. Fix xq S suppA(io)- Let e C^{Bp{x{))) with (j){xo) ^ 0. 
Then if I is sufficiently short, we find 

0dA(i) 7^0 
'n 

for t ^ I. Hence Bp{xo) intersects suppA(i). It follows that ^{t) > if (to), and by 
the choice of to, this implies equality. D 

We now fix the interval / from Lemma [6.21 and we set Q — #(to)- We may 
relabel the points bp{t) such that 

Q 

Ht) ^E'^pW^^w- 

p=i 

The continuity of X{t) and the fact that a-p{t) G 47rZ\{0} then imply that bp{t) is 
continuous in / (if labeled appropriately) for p = 1, . . . , Q. 

Lemma 6.3. For almost all t £ I and all (f) G C^{ilr{a{ti))), 

p 

J2 V^V0(6p(i)) : Bp{t) - 0. 
P=Q+1 

Proof. It suffices to prove this for a with bp{t) ^ supp (j) for p ~ 1, . . . ,Q. But in 
this case, the function 

Q 

ih^^f7p(i)0(&p(i)) 
p=i 

is constant in a neighborhood of t. The claim then follow from (|20p . D 

Lemma 6.4. /n (ii,i2); 
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Proof. If we insert functions (f> in (pn| with V^^(6p(t)) — for p = 1, . . . ,Q at a 
time t G I, then by density arguments we obtain the required identity in /. But 
then the curves bi{t), . . . , bQ{t) arc constant in I, and we can extend the arguments 
beyond the end points of / (unless they coincide with ti or ^2)- D 

The remaining arguments are similar to [52] . Now consider (IT5|) again and insert 
(j),ip € C^i^) with V^(/) = V-0 and V^V- = in 17\f7r(a(ii)). Then we obtain 

N 

TT^ (aoV'(an(^2)) +4g„(t2)0(a„(f2)) -aoV'(an(^i)) - 4g„(ti)(?!)(a„(ti))) 

n=l 

= f[ f V-^Vf/) : rfAi*(^) + 51 V-LV(/)(6p(i)) : Bp(i) j dt 

for certain functions qi, . . . ,qN : [ti,t2] ^> 2 + '^, using Lemma 3 in [22] again. 
If (f> — ^j — in in n\Qrio,{ti)), then we have the identity 

Jti Jn 
again. It follows that 

P rt2 



^ / V^ V0(5p(t)) : Bp{t) dt = 



for all (f), ijj with this property. But then the last identity must be true for all 
<j>,tp & C^{il.). Hence we have in fact 



N 



TT ^ {aQip{an{t2)) + 4g„(t2)0(a„(i2)) - aotlj{an(ti)) - 4g„(ii)0(a„(ii))) 

n=l 

= / / V^\7(j) : dfi4t)dt. 
Jti Jn 

for all 4>,ip & C^{fl) with V-^cj) — Vtp in a neighborhood of the vortices. From 
this we see that each (?„ is locally constant away from E. Using p^ . we derive the 
motion law 

d 

(ao + 4g„i)7rd„ = -- — W{a, d), n = l,...,N. 
dan 



Thus we have proved Theorem 1 1.1 



7. Comments on the complex valued case 

As remarked in the introduction, a similar analysis can be performed on the 
complex Ginzburg-Landau equation of mixed-type dynamics, 

(21) {a,+t)^ = Au, + -u,{l-\u,\^) 

on XI C M^ with ae log ^ — >■ ao- The dynamics of vortices were established rigorously 
by Miot [55] on K^ and by the authors [21] in bounded domains. In both papers the 
authors make use of compactness results for the Ginzburg-Landau energy egi(uc) 



d , , 




du 


^ f du 


9^^^(") = 


-OLe 


dt 


+ div Vu, -— - 
\ dt 
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and the Jacobian J{u^), the associated conservation laws 

(22) |0(H^-l))=divjK)-«.(z.,|^) 

( Ou 

(23) Y*^*^") " curl div (Vu (g) Vu) + a^ I -^, Vu 

(24) 

and well-preparedness of the initial data to pin vortices to the ODE. 

In order to prove Theorem 11.21 we essentially follow the arguments from the 
proof of Theorem ll.il however, there is one important improvement. Solutions of 
equation (I2ip remain smooth for all times, and thus there is no bubbling and no 
discontinuity in the coefficients of the ODE. Furthermore, due to the compactness 
results of Jerrard-Soner 'T^ and continuity of the Jacobian, one finds that J{Uf) — >■ 
J{u^). In the LLG case, on the other hand, we cannot prove uj[mf:) — > a;(m,), as 
there may be additional Dirac masses induced by bubbling. Consequently, equation 
(PT|) permits a few shortcuts in the proof. Otherwise, the arguments remain the 
same and we do not repeat them. 

The procedure follows along the lines of the proof of Theorem 11.11 First we 
establish strong convergence of the stress energy tensor ,^'} (g) ''| " | to jiu<t)(g)j{u^) 
outside the vortex cores and a finite number of points. The proof is essentially the 
same as for Lemma [Ql and uses identity ([2^ . Combining the differential identities 
(P^ and (IMl) , in a similar way as in the proof of Theorem 11.11 we show that 
the defect measure does not affect the motion of the vortices, and the motion law 
follows. 
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